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Abstract We analyze the 26 November 2005 solar radio event observed in- 
tcrfcromctrically at frequencies of 244 and 611 MHz by the Giant Metrewave 
Radio Telescope (GMRT) in Pune, India. These observations are used to make 
interferometric maps of the event at both frequencies with the time cadence of 1 s 
from 06:50 to 07:12 UT. These maps reveal several radio sources. The light curves 
of these sources show that only two sources at 244 MHz and 611 MHz are well 
correlated in time. The EUV flare is more localized with flare loops located rather 
away from the radio sources. Using the SoHO/MDI observations and potential 
magnetic field extrapolation we demonstrate that both the correlated sources 
are located in the fan structure of magnetic field lines starting from a coronal 
magnetic null point. Wavelet analysis of the light curves of the radio sources 
detects tadpoles with periods in the range P = 10-83 s. These wavelet tadpoles 
indicate the presence of fast magnetoacoustic waves that propagate in the fan 
structure of the coronal magnetic null point. We estimate the plasma parameters 
in the studied radio sources and find them consistent with the presented scenario 
involving the coronal magnetic null point. 
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1. Introduction 



It is commonly believed that the primary energy release regions in solar flares 



arc located in the low corona (Aschwanden, 20051. Radio spectral and imaging 



observations in the dccimctric and metric wavelength ranges in combination 
with magnetic field extrapolations are considered to be very promising tools for 
a study of these processes. Unfortunately, imaging observations of solar flares in 
the decimetric range arc very rare. In the metric range, such observations are 



commonly done by the Nangay radiohcliograph (Kerdraon and Delouis, 19971. 

To successfully combine radio observations and magnetic extrapolations, the 
radio observations have to include sufficiently precise positional information. 
Such a combination was done e.g. bv ITrottet et al.\ (|2006p . where a detailed 
analysis of radio spectral and imaging observations in the 10-4500 MHz range 
was presented for the 5 November 1998 flare. ISubramanian et~al] (j2007|) have 
studied a post-fiarc source imaged at 1060 MHz to calculate the power budget 
for the efficiency of the plasma emission mechanism in a post-flare decimet- 
ric continuum source. |Aurass, Rausche, and Mann| (I2007P have analyzed the 
topology of the potential coronal magnetic flcld near the source site of the metcr- 
decimetcr radio continuum to flnd that this radio source occurs near the contact 
of three separatrixes between magnetic flux cells. lAurass et al. \ (|201ip have ex- 
amined meter-decimeter dynamic radio spectra and imaging with longitudinal 
magnetic field magnetograms to describe meter- wave sources. [Chen et a/.l(|201ip 
have used an interferometric dm-radio observation and nonlinear force-free field 
extrapolation to explore the zebra pattern source in relation to the magnetic 
field configuration. Zuccare llo et al.\ (|201ip investigated the morphological and 
magnetic evolution of an active region before and during an X3.8 long duration 
event. They found that coronal magnetic null points played an important role 
in this flare. 

A comprehensive review of magnetic fields and magnetic reconnection theory 
as well as observational findings was provided bv lAschwandenI (|2005p . The topo- 



logical methods for the analysis of magnetic fields were reviewed in Longcope] 



(|2005p and the 3D null point reconnection regimes in lPriest and PontinI (|2009p . 
[McLaughlin, Hood, and De Moortel] (j201ip have presented a review of the the- 
oretical studies of the MHD waves in the vicinity of the magnetic null points. 
Furthermore, [Afanasyev and Uralov (j2012p have studied analytically the prop- 



agation of a fast-mode magnetohydrodynamic wave near a 2D magnetic null 
point. Using the nonlinear geometrical acoustic method they have found complex 
behavior of these waves in the vicinity of this point. In spite of the wealth of the- 
oretical work presented in these papers, the authors concluded that there is still 
no clear observational evidence for the presence of MHD waves near null points of 
the magnetic field. We note that this is also in spite the fact that MHD waves are 
commonly observed in the solar corona Pe Moortel, Ireland, and Walsh, 20001 



Kliem et al, 2002[ [Harrison, Hood, and Pike, 2002l[Tomczyk et ai, 2"007|[Ofman and Wang, 2008H 



De Moortel, 2009 Marsh, Walsh, and Plunkett, 2009 Marsh, De Moortel, and Walsh, 20TT1).| 



Since both the standing and propagating magnetoacoustic waves modulate the 
plasma density and the magnetic field in the radio source ( [Aschwanden, 2005^ , 
some modulation of the radio emission by both these waves can be expected. 
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Figure 1. Solar maps with main sources associated with the 26 November 2005 radio event 
observed at 7:01:30 UT by the GMRT instrument. Regions U1-U3 and D1-D3 were identified 
as the main radio sources at 244 and 611 MHz, respectively. The magenta circle in the left 
panel has the diameter of 32 arc min and indicates an approximate position and size of the 
visible solar disk. Disk centers are shown by a cross. The cross has dimensions of 400" in both 
directions to indicate the scale in the maps. Synthesized beam dimensions giving the error in 
GMRT positions are represented by the small ovals shown on the bottom right corners. 



Roberts, Edwin, and Benz (|19831 I1984p studied impulsively generated fast 



magnetoacoustic waves trapped in a structure with enhanced density {e.g. loop). 
They showed that these propagating waves exhibit both periodic and quasi- 
periodic phases. INakariakov et aL\ (|2004p numerically modelled impulsively gen- 
erated fast magnetoacoustic wave trains and showed that the quasi-periodicity 
is a consequence of the dispersion of the guided modes. Using wavelet analy- 
sis, these authors found that typical wavelet spectrum of such fast magnetoa- 
coustic wave trains is a tadpole consisting of a broadband head preceded by 
a narrowband tail. 

The tadpoles as characteristic wavelet signatures of fast magnetoacoustic wave 
trains were observed in solar eclipse data ( Katsiyannis et ai, 2003 ) as well as in 
radio spectra of decimetric gyrosynchrotron emission ( [Meszarosova et ai, 2009a| ),| 
and also in decimetric plasma emission ( Meszarosova et ai, 2009b I. While the 
tadpoles in the gyrosynchrotron emission was detected simultaneously at all ra- 
dio frequencies, the tadpoles in the plasma emission drifted towards low frequen- 
cies. This type of "drifting tadpoles" was studied in details in |Meszarosova, Karlicky, and Rybak[ 
(|201ip in the radio dynamical spectrum with fibers. 

The observed parameters of fast magnetoacoustic waves reflect properties of 
the plasma in the waveguides where these waves are propagating. Therefore, 
one could use observed waves and their wavelet tadpoles as a potentially useful 
diagnostic tool ( [Jelfnek and Karlicky, 2009 Jelfnek and Karlicky, 20l"0| for de- 
termining physical conditions in these waveguides, {e.g. loops or current sheets). 
[Karlicky, Jelfnek, and Meszarosovd (j201ip compared parameters of wavelet tad- 
poles detected in the radio dynamical spectra with narrowband spikes to those 
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Table 1. Time evolution of the individual GMRT 
sources. 



Source Start time Time of max End time 



[UT] [UT] [UT] 



Ul 


6:58:04 


7:01:30 


7:04:03 


U2 


6:57:59 


7:01:30 




U3 


6:58:20 


7:06:45 




Dl 


6:58:05 


7:03:52 


7:05:11 


D2 


6:58:54 


6:59:33 


7:04:00 


D3 


6:58:13 


6:58:43 


7:01:59 



computed in the model with the Harris current sheet. Based on this comparison 
the authors proposed that the spikes are generated by driven coalescence and 
fragmentation processes in turbulent reconncction outflows. We note here that, in 
general, flare current sheets can be formed not only near magnetic null points, 
but also e.g. between interacting magnetic loops. |Jehnek and Karlicky| (12012^ 
numerically studied impulsively generated magnetoacoustic waves for the Harris 
current sheet and a density slab. In both cases they find that wave trains were 
generated and propagated in a similar way for similar geometrical and plasma 
parameters. 

In this paper, we analyze a rare decimetric imaging observation of the 26 Novem-| 
ber 2005 solar flare made by the Giant Metrewave Radio Telescope (GMRT). 
Combing the results of this analysis with the magnetic field extrapolation (Sec- 
tion[2]), we present a scenario of this radio event. For the first time, we detected 
the magnetoacoustic waves in the radio sources (Section located in the fan 
of magnetic field lines connected with a coronal null point. The basic plasma 
parameters in the radio sources are estimated and the results are discussed 
(Section El). 



2. Observations and data analysis 

The B8.9 solar fiare occurred on 26 November 2005 in the active regions NOAA 
AR 10824 and 10825 located near the disk center. The fiare lasted from approx- 
imately 06:31 to 07:49 UT with the GOES maximum at 07:05 UT. 

The radio counterpart of this flare was a 22 minutes long radio event lasting 
from the 06:50 to 07:12 UT recorded by the Giant Metrewave Radio Telescope 
(GMRT) in Pune, India. The Michelson Doppler Imager (MDI, IScherrer et al.\ 
I1995P onboard the SoHO spacecraft ( [Domingo, Fleck, and Poland, 1995[ ) per- 
formed routine full-disk observations with a 96 minute cadence. The magne- 
togram nearest to the radio event was observed at 06:24 UT. The EUV counter- 
parts of this flare were observed by the Extreme Ultraviolet Imaging Telescope 
dPelaboudiniere et al., 1995] ) onboard the SoHO spacecraft. 
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Figure 2. Selected contour maps showing the time evolution of the emission sources at 
244 MHz (upper panel) and 611 MHz (bottom panel). The times (UT) are in blue. Synthesized 
beam dimensions representing the error in GMRT positions are shown as small green circles 
on the bottom left corners of maps at 6:57 UT. Upper panel: Positions of the sources U1-U3 
(in red) are indicated on the map at 6:58:30 UT. Bottom panel: Positions of the sources D1-D3 
(in red) are indicated on the map at 6:59 UT. 



2.1. GMRT observations and data analysis 

On 26 November 2005, the GMRT observed the Sun at two frequencies, 244 and 
611 MHz. The GMRT instrument ( |Swarup et al., 1991|[Ananthakrishnan and Rao, 2002H 
|Mercier et ai, 2006[ ) is a radio interferometer consisting of 30 fully steerable indi- 
vidual radio telescopes. Each telescope has a parabolic antenna with a diameter 
of 45 m and 16 of these individual antennas are arranged in an Y shape array 
with each arm having length of 14 km from the array centre. The remaining 14 
telescopes are located in the central area of 1 km^ . The interferometer operates 
at wavelengths longer than 21 cm, with six frequency bands centered on the 38, 
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Figure 3. Light curves of the 26 November 2005 radio event obtained from the GMRT maps. 
Upper panel: Radio fluxes at 244 MHz corresponding to the sources Ul (blue), U2 (dashed 
green) and U3 (dotted red). Bottom panel: Radio fluxes at 611 MHz corresponding to the 
sources Dl (blue), D2 (dashed green) and D3 (dotted red). 

153, 233, 327, 610, and 1420 MHz frequencies. The maximum resolution depends 
on the configuration, and varies between 2 and 60 arc sees. 

The observed interferometric data at 244 and 611 MHz were Fourier trans- 
formed to generate a series of 1320 snapshot images of the Sun at 1 second 
time-cadence, from 06:50 to 07:12 UT. The images were cleaned using the algo- 
rithm developed by I Schwab I (|1984p and rotated to correct for the solar North. 
The synthesized beam dimensions giving the GMRT positional error are are 
77.7 X 50.8 and 17.7 x 13.4arcsec at 244 and 611 MHz, respectively. 

An example of the observations showing the main radio sources is shown in 
Figure[TJ On this figure, six sources can be identified, three for each frequency. At 
the frequency of 244 MHz, the sources are outlined in the left panel and denoted 
as Ul, U2, and U3. The main sources at the frequency 611 MHz are shown in the 
right panel of Figure [T] and labelled as Dl, D2, and D3. The evolution of these 
sources during the radio event is shown in Figure [5] Positions of the individual 
radio sources U1-U3 and D1-D3 arc shown in red. Notable is the merging of 
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Figure 4. Superposition of the radio sources at 244 MHz (orange contours) and 611 MHz 
(blue contours) observed at 06:59 UT on top of the photospheric magnetic map obtained with 
the SoHO/MDI at 6:24 UT. Individual radio sources are labelled. The MDI magnetogram is 
saturated to ± 10"^ G. Heliographic coordinates are shown for orientation. Solar North is up. 

the sources Ul and U2 at around 6:59 UT. The source U3 remains far from the 
others (Ul and U2). 

We constructed hght curves of the individual radio sources by enclosing these 
sources in rectangular regions on the GMRT maps. The light curves are shown 
in Figure 131 The radio fluxes at 244 MHz corresponding to the sources Ul, U2, 
and U3 are shown in the upper panel as the blue, dashed green and dotted red 
lines, respectively. The radio fluxes at 611 MHz frequency corresponding to the 
sources Dl (blue), D2 (dashed green) and D3 (dotted red) arc present in the 
bottom panel. The calibration of the individual fluxes is relative only and is 
expressed in arbitrary units (a.u.). 

Analysis of these time scries of radio fluxes shows different time evolution 
profiles for the different radio sources. The sources Ul, Dl, and D3 show about 
5 minutes durations with a well defined main peak. On the other hand, the 
sources U2, U3, and D2 show profiles with a gradual rise and several peaks on 
top of a roughly constant background. The temporal properties of these sources 
are presented in Table 1. The start and end times of each source were determined 
as the times when the fluxes were above or below half the average flux of the 
whole burst profile of the source, respectively. 
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Figure 5. Superposition of the radio sources (same as in Figure |4]l on top of the SoHO/EIT 
195 A observations taken at 7:48 UT. The BIT 195 A observations are saturated to 2x10"^ 
DNs~^ px~^ and shown in logarithmic intensity scale for better visibility of fainter emitting 
loop systems. The flare arcade footpoints observed using the filter 304 A at 07:19 UT arc showrn 
as white contours corresponding to 10"^ DNs~^ px"'^. 



The cross-correlation coefficients between pairs of all GMRT U and D sources 
show that there is only one pair with correlation higher than 75%. This pair is 
the sources Ul and Dl with the cross-correlation coefficient of 81%. This degree 
of correlation indicates that these radio sources can have a common origin. 

The maximum of the cross-correlation coefficient is flat if the time lag of Ul 
with respect to Dl ranges 0-50 s. It indicates that first parts of radio fluxes Dl 
and Ul are correlated by some fast moving agents (possibly beams) and the 
second parts by slowly moving agents (possibly waves). 

2.2. Photosphcric magnetic field and the EUV flare 

The GMRT interferometric observations offer the advantage of direct spatial 
comparison with the photosphcric magnetic field and the EUV flare morphology, 
as observed by the SoHO/MDI and SoHO/EIT instruments, respectively. The 
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Table 2. Basic parameters derived for radio bursts at GMRT frequencies. 



Frequency Source 
[MHz] 



Fundamental 
frequency altitude 
[Mm] 



Plasma density 
[cm^''] 



First harmonic 
altitude 

[Mm] 



244 
611 



U1-U3 
D1-D3 



48 
22 



7.4 X 10** 
4.6 X 10^ 



40 



comparison of the radio sources loci with the JVIDI magnetograni is presented in 
Figure m The magnetogram was observed at 06:24 UT and rotated to the time 
06:59 corresponding to the radio observations using the SolarSoftware routine 
drot_map.pro. The time of the radio observations, 06:59 UT, is chosen because it 
corresponds to the times of the first radio peaks in the sources Ul and Dl. 

The radio sources are located in a quadrupolar magnetic configuration consist- 
ing of two active regions, NOAA 10824 and 10825 (FigureS]). Both active regions 
arc bipolar with a /3-configuration. The AR 10824 is located approximately at the 
heliographic latitude of ~ —13° and contains a wcU-dcvelopcd, leading negative- 
polarity sunspot. Other polarities consist of plage regions or pores, which is the 
case also for the AR 10825 located at latitudes of ~ —7°. A notable feature is 
that the radio sources arc not located on top of the main magnetic polarities, 
but, in general, they overlie weak-field regions. 

The SoHO/EIT instrument was performing full-disc observations in the 195 A fil-| 
ter with a cadence of approximately 12 minutes. Observations in other filters, i.e. 
171, 284 and 304 A were performed at 07:00, 07:06 and 07:19 UT, respectively. 
Due to poor pointing information of the EIT instrument, the EIT 304 A ob- 
servations were coaligned manually with the JVIDI observations by matching 
the EIT 304 A brightenings with small magnetic polarities observed by JMDI, 
which show good spatial correlation ( [Ravindra and Venkatakrishnan, 2003[ ). We 
estimate the error of this manual coalignment to be w 5 arc sec. 

In the EIT 304 A observations, three fiare arcade footpoints are discernible, 
shown in Figure [5] as white contours corresponding to the observed intensity 
of 10^ DNs^^px"^ AU three footpoints are located weU within the AR 10824, 
with one of them in the positive polarity and the other two in the negative 
polarities North of the sunspot. EIT 195 A observations show cooling system of 
flare loops connecting these three footpoints. The flare loops are well- visible at 
195 A at 07:48 UT and are shown in Figure [S] The global magnetic configuration 
of AR 10824 is that of a sigmoid with large shear. The magnetic configuration 
of the neighboring AR 10825 is near-potential. 

Comparison of the loci of the radio sources with the EUV fiare morphology 
(Figure [S]) shows that the radio sources have no direct spatial correspondence 
with the EUV fiare loops or their footpoints. The source D3 is an exception, 
since it overlies a portion of the flare loops. The sources Ul, U3, Dl and D2 are 
located in the area of weak EUV emission. 
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3. Magnetic structure of active regions NOAA 10824 and 10825 

3.1. Magnetic field extrapolation and the altitude of the radio emission 

To investigate the relationship between the radio sources and the structure of 
the magnetic field of active regions NOAA 10824 and 10825, we performed an 
extrapolation of the SoHO/MDI magnetogram (Figure H]) observed at 06:24 UT 
prior to the flare and associated radio events. The extrapolation was carried 
out in linear force-free approximation, where the magnetic field B given by the 
solution of the equation 

VxB^aB, (1) 

for a = const. The solution is subject to the boundary condition Bz{x, y, z ~ 0) 
given by the observed magnetogram, where < a; < i^; and < y < Ly. The 
constant a is subject to the condition a < amax = 27r/max(L2;, Ly), otherwise 
the magnetic field is non-physical. We utilized the Fourier transform method 
developed by lAlissandrakisI (|198ip and |Gary| (|1989p . This method allows for 
extrapolation of the part of the observed magnetogram in a carthesian geometry. 
The computational box is shown in Figures [6] and [T] 

We calculated a range of linear force-free models with various values of a. 
However, the flare loops are poorly approximated with a = const., even with 
large values of a close to amax- The reason for this probably is the presence 
of differential shear within the active region ( [Schmieder et al., 1996[ ). Moreover, 
using large values of a leads to poor flt to the observed shape of the coronal loops 
in the AR 10825, which is close to the potential state (a = 0). Therefore, we chose 
to extrapolate in the potential approximation. We also note that the potential 
approximation does usually a good job in capturing the topological structure of 
the active region (though not at sigmoid locations, e. a. iSchmieder and Aulanierl 

mm . 

To compare the 3D magnetic held geometry with the loci of the radio sources, 
observed in a 2D plane of the sky, the approximate altitude at which the radio 
emission originates must be determined. To do that, we consider the parameters 
of the radio bursts (at GMRT frequencies) and the solar atmosphere density 
model ( |Aschwanden, 2002] ) for the radio plasma emission at fundamental fre- 
quency and the first harmonic (see also Section 13. 2p . We use Aschwanden's 
density model because it was derived from radio observations. The basic pa- 
rameters of the radio sources are derived and summarized in Table [5] where the 
plasma density values belong to the fundamental frequency altitude. 

3.2. Magnetic topology and its relation to radio sources 

The extrapolated magnetic field contains a pair of coronal magnetic null points 
between the two active regions. The negative and positive coronal null points are 
denoted by red (SE) and blue (NW) arrows in Figure HI On this figure, the dark- 
red and dark-blue contours correspond to positive and negative photospheric 
polarities, respectively, with B^ ~ ±50 and 500 G. Since the magnetic field of 
AR 10825 is weaker, the null points are located closer to this active region than 
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Figure 6. Plane-of-the-sky projection of the extrapolated magnetic field. Red (SE) and blue 
(NW) arrows show the locations of the positive and negative coronal magnetic null points, 
respectively. The separator lying at the intersection of the respective fans is plotted as thick 
red line. The magnetic field lines are colored according to the local magnetic field, in the 
range of logio(B/T) S (— 1.0, 3.5). The thin, dark-red and dark-blue contours denote posi- 
tive and negative photospheric polarities, respectively. The two contour levels correspond 
to Bz = 0) = ± 50 and 500 G. Shades of gray show the photospheric intersections of 
the quasi-separatrix layers. The contours of the GMRT sources are located at the altitudes 
corresponding to the fundamental frequency (Table [JJl and are shown as thick orange (Ul-3) 
and blue (Dl-S) contours, respectively. 



to the larger AR 10824. The existence of these null points is a consequence of 
the multipolar structure of the magnetic field {e.g. |Longcope|20051 and references 
therein). The magnetic structure in the vicinity of these null points is shown us- 
ing several field lines regularly distributed near the null points. The fan surfaces 
of these null points intersect, forming a separator ( |Baum and Bratenahl, 1980[ 
|Longcope, 2005[ ). The separator is found manually by trial-and-error and is 
shown by a thick red line. 
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Figure 7. Same as in Figure[6] but in a 3D projection to show tlie altitudes of the individual 
field lines as well as individual radio sources. Only the field lines lying in the vicinity of the 
negative null point or passing through the source D2 are shown for clarity. 



Both the spine and the fan surface of the negative null point are closed within 
the computational domain, while the fan surface of the positive null point is open. 
This is not typical of coronal null points ( [Pariat, Antiochos, and DeVore, 2009{ 



Del Zanna et al., 2011 ). Both fan surfaces form a part of the quasi-separatrix lay- 



ers ( Priest and Demoulin, 1995{|Demoulin et al, 1996[|Demoulin et al, 1997| ) as-| 



sociated with the active regions. The intersections of these quasi-separatrix layers 
with the photosphere are shown in Figures |6] and [7] by shades of gray corre- 
sponding to the unsigned norm of the field line footpoint mapping, defined 
as 




Ni^^y)- ^ + -J- + -J- + -J- > (2) 
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where the starting footpoint of a field hnc at z = is denoted as {x, y) and 
its opposite footpoint as {X,Y). The unsigned norm is computed irrespective 
of the direction of the footpoint mapping. We note that the correct definition 
of quasi-separatrix layers is through the squashing factor Q or the expansion- 
contraction factor K ( |Titov, Hornig, and Demoulin, 2002 ), which are invariants 
with respect to the footpoint mapping direction. Here, we use the unsigned norm, 
since it is much less sensitive to the errors arising because of the measurement 
noise and numerical derivation due to finite magnetogram resolution, and at the 
same time gives a good indication of the expansion-contraction factor K . 

The sources U and D are plotted as thick blue and thick orange contours, 
respectively. These sources are plotted in height corresponding to the altitudes 
given by the fundamental frequencies (Table [2]). It is clearly seen that both Ul 
and Dl are lying along the fan of the negative null point, with the field lines 
passing through the Dl being separated by the separator (Figures [5] and [7]). 
Therefore, both the Ul and Dl lie along a common magnetic structure, provid- 
ing straightforward explanation for the observed correlation between these two 
sources (Section [^?T|) . The separator also passes through Ul, as does some other 
fan field lines rooted in AR 10825. 

We note that the fan of the negative null-point covers only a part of the source 
Ul. The rest of the Ul is associated with closed field lines rooted in the positive 
polarity of the AR 10825, in direct vicinity of the footpoints of the fan field 
lines. One example of such field lines is plotted in Figures [S] and [T] That there 
are structures other than the fan of the negative null-point passing through the 
source Ul offers explanation for the fact that the correlation coefficient between 
Ul and Dl is smaller than one. 

The sources D2 and U3 are located on the field lines that are open within the 
computational box. The footpoints of these field lines arc located in the AR 10824 
in the vicinity of the photosphcric quasi-separatrices with highest A'^ separating 
the open and closed magnetic flux (Figure (5]). Overall, we conclude that the field 
lines passing through the sources Dl, D2, Ul and U3 arc rooted in the vicinity 
of the flare arcade footpoints (compare Figures [5] and (B]) , while the source D3 
appears to be connected directly to flare loops. These observations point out 
the fact that even a localized, weak flare can cause a widespread perturbation 
involving a large portion of the surrounding magnetic field and its topological 
features. 

We stress that the goodness of the match between the extrapolated magnetic 
field and the radio sources is subject to the approximations used. E.g. presence 
of electric current in the model of the magnetic field would alter the geometry of 
the magnetic field. In particular, using a linear force-free approximation would 
mostly affect the longest field lines passing through the source Ul. Additionally, 
the density model of ( |Aschwanden, 2002] ) provides only a single altitude for the 
emission of each radio frequency, while in reality the density profile and thus the 
altitude of emission can be different for each field line. However, the available 
longitudinal MDI magnetogram and the KIT data do not provide necessary 
constraints (Section ^ for using more elaborate models of the magnetic field or 
density throughout the observed radio-emitting atmosphere. 
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Table 3. Characteristics periods P for 
wavelet power spectra showing a tadpole 
pattern. 





Source 


Frequency 
[MHz] 


Period P 


Ul 
Dl 


244 
611 


83, 25, 22, 11, 10 

76, 22, 13 





Figure 8. Examples of the wavelet power spectra showing tadpole patterns with characteris- 
tics periods P (bottom panels) in comparison with the GMRT radio fluxes (upper panels) at 
244 MHz (source Ul) and 611 MHz (source Dl) frequencies. The lighter area shows a greater 
power in the power spectrum and it is bound by the solid contour of the confidence level >95%. 
The hatched regions belong to the cone of influence with edge effects due to finite-length time 
series. The arrows (1) and (2) correspond to 7:01:30 and 7:03:52 UT, respectively. 



We also note that the magnetic connectivity between Ul and Dl exists only 
if these sources are assumed to emit at the fundamental frequency (Tabic [5]) . 
If the origin of the radio emission would be the second harmonic, the radio 
emission should occur at altitudes which are higher than the vertical extent of 
the negative null point fan. This point is discussed further in Section [51 



4. Wavelet analysis of the GMRT time series 



To investigate the nature of the disturbance giving rise to the radio sources, an 
analysis of possible periodicity of radio time series at all frequencies, detected by 
the GMRT instrument, was made using the Morlet wavelet transform method 
( Torrence and Compo, 1998 1. We focused on the well correlated GMRT sources 
Ul and Dl. 

The wavelet power spectra at all frequencies of the GMRT instrument show 
the tadpole patterns ([Nakariak ov et al., 2004[ ) with characteristics periods P 
([Meszarosova et aZ.i [Migau2009b| in the range 10-83 s. 

These periods were detected in the entire studied time interval 06:50-07:12 UT 
and arc listed in Tablc[3] The periods of about 10 s were found at all frequencies. 
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Dl 



Figure 9. Scheme of the suggested scenario for the well-correlated sources Ul and Dl lying 
along the fan of the coronal null point. Magnetoacoustic waves move along magnetic field lines 
(trajectories (1) and (2)). The separator is located between these trajectories. 

Figure [5] shows two examples of the wavelet power spectra of the well correlated 
radio sources Ul and Dl showing tadpole patterns (bottom panels) with the 
longest characteristic periods P listed in Table [S] In this figure, the arrows (1) 
and (2) are used to denote the times of 7:01:30 and 7:03:52 UT, respectively. 
These times correspond to the main peaks of the flux shown in the upper panels 
and the occurring tadpole head maximum shown in bottom panels. It has been 
proposed that the tadpoles indicate arrival of the magnetoacoustic wave trains 
at the radio source (jNakariakov et al.[ 120041 IMeszarosova et al.[ I2009al I2009bp . 

Although the Ul and Dl radio fluxes are well correlated (81%) the sources 
Ul and Dl are not the same in all details. Different appearance of the tad- 
poles in Figure H] reflects different plasma parameters in the sources Ul and Dl 
( [Jelmek and Karlicky, 2012] ). This is expected, since these tadpoles have differ- 
ent characteristic period P, which means that they are propagating in different 
waveguides (see Section [S]). 



5. Discussion and Conclusions 

We studied solar intcrferomctric maps (Figurcs[T]and[5]) and time scries of the six 
individual radio sources, U1-U3 at 244 MHz and D1-D3 at 611 MHz (Figure[3l), 
observed by the GMRT instrument with 1 s time resolution during the B8.9 flare 
on 26 November 2005. We determined that only sources Ul and Dl arc well 
correlated with the cross-correlation coefficient of 81%. 

Generally, the positions of the radio sources can be influenced by refractive 
effects. They play a role in the solar atmosphere as well as in the Earth iono- 
sphere. As concerns to the effects in the solar atmosphere, our low frequency 
244 MHz source Ul (which is the main feature under study together with the 
source Dl) is nearly in the disk center, see Figure [T] or |4l In the disk center the 
refractive effects come to zero. Similarly, the Dl source at 611 MHz is close to 
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the disk center. We note that the phase cahbration process at GMRT corrects 
the ionospheric refraction. 

Comparison of the GMRT interferometric maps with the EUV observations 
of the flare by the SoHO /EIT instrument showed that only the source D3 can be 
connected to the system of flare loops. All other radio sources are located well 
away from the X-ray and EUV flare (c/. |Bcnz, Sa int- Hilaire, and Vilmer [2002p . 
Potential extrapolation of the observed SoHO/MDI magnetogram in combina- 
tion with the solar coronal density model of I AschwandenI (|2002p allowed us to 
find the connection of these other radio sources to the magnetic field of the 
observed active regions and their magnetic topology. The sources Ul and Dl 
were found to lie along the fan of a negative coronal null point. That the sources 
Ul and Dl lie along a common magnetic structure offers explanation for their 
observed correlation. The sources D2 and U3 lie at open field lines anchored near 
the quasi-separatrices in the vicinity of the flare arcade footpoints. 

The connection between the radio sources and the magnetic field, as found 
in Section 13.21 is valid only if the radio emission originates at the fundamental 
frequency. This is because the altitudes corresponding to the first harmonic are 
too high (Table [5]) . We note that the emission at the first harmonic is usually 
considered to be stronger than at the fundamental frequency, due to the strong 
absorption for the fundamental emission. However, there are cases {e.g. enhanced 
plasma turbulence in the radio sources) which reduces this absorption and thus 
the emission at the fundamental frequency can be stronger than at the first 
harmonic. 

Based on the results above, we conclude that even the observed localized 
fiare is able to cause a widespread perturbation involving a large portion of the 
topological structure of the two active regions. This perturbation can be under- 
stood in the following terms: As the flare progresses, the flare arcade footpoints 
move away from each other. Since these footpoints correspond to the intersec- 
tions of the quasi-separatrix layers with the photosphere ( [Demoulin et aL, 1996[ 



Demoulin et at, 1997), the increase of their distance perturbs the surround- 



ing magnetic field and excite waves. If the surrounding, perturbed magnetic 
field contains null points, these can collapse and form a current sheet, thereby 
commencing reconnection in regions not directly involved in the original fiare. 
This can then be also the mechanism for sympathetic fiares and eruptions 
dMoon et aZ.7^ 02: Kha n and Aurass, 2006[|Jiang et al, 2008[[Torok et al, 2011H 
|Shen, Liu, aSdSu, 2012| ). 

For the first time, wavelet tadpole patterns with the characteristic peri- 
ods of 10-83 s (Table [3]) were found at metric radio frequencies. We have in- 
terpreted them in accordance with the works of INakariakov et al.\ (|2004p and 
IMeszarosova et aD (j2009al I2009bp as signatures of the fast magnetoacoustic 
waves propagating from their initiation site to studied radio sources Ul and 
Dl. 

The mechanism for initiation of the inferred magnetoacoustic waves present in 
the radio sources can thus be both the observed EUV flare and the flare-induced 
collapse of the null point located in the surrounding magnetic field. These waves 
propagated towards radio sources along magnetic field lines. They arrived at 
the radio sources (sec the wavelet tadpoles in Figure and modulated the radio 
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emission there. We expect that the magnetoacoustic waves, through their density 
and magnetic field variations, modulate the growth rates of instabilities (like the 
bump-on-tail or loss-cone instabilities), which generate plasma- as well as the 
observed radio waves. 

The schematic scenario for the well-correlated sources Ul and Dl lying along 
the fan of the coronal null point is shown in Figure IHl Magnetoacoustic wave 
trains move along magnetic field lines (trajectories 1 and 2). The separator 
is located between these trajectories (i.e. between two groups of magnetic field 
lines (c/. Figure [5]). The distances between the null point and radio sources along 
the magnetic field lines are obtained directly from the extrapolation and allow 
the determination of the velocities of the magnetoacoustic waves. The averaged 
distances between the radio source and the null point, are about 53 and 103 Mm 
for the Ul and Dl, respectively. We considered the time intervals t2 — tl and 
tS — tl where tl is the time of triggering of the magnetoacoustic wave trains. We 
assume that these trains propagating to the source Ul as well as to the source 
Dl were generated simultaneously at the beginning of the radio event, i.e. at the 
start time of Ul and Dl sources (Table [Ij. The times t2 and ^3 correspond to 
the times of tadpole head maxima (arrows (1) and (2), Figure respectively. 
Thus, the mean velocities of the magnetoacoustic waves are 260 and 300 km s^^ 
for Ul and Dl sources, respectively. 

Now, let us compare these mean velocities with the Alfvcn velocities at the 
radio sources at Ul and Dl. Taking plasma densities from Table and the 
magnetic field from the extrapolation (_B = 3G at Ul and 13 G at Dl) the 
Alfven velocity at these sources are va = 220 and 390 km s"-'^, respectively. These 
velocities are in agreement with previous results, considering that the Alfven 
velocities change along the trajectory of these magnetoacoustic waves, starting 
from the coronal null point. 

Thus, knowing the Alfven velocities in the radio sources, we can estimate also 
the width w of the structure through which the magnetoacoustic wave trains 
are guided. Namely, the period of this magnetoacoustic wave can be estimated 
as P ~ w/vA (Nakariakov et al. , 2004). Thus, the widths w of the struc- 
tures guiding the magnetoacoustic waves (modulating the radio emission) are 
18Mm (83s x 220 kms'^) and 30 Mm (76s x 390 kms^^) in the Ul and Dl 
radio sources, respectively. This rough estimation confirms that the extent of 
the structure guiding magnetoacoustic waves is within the width of the fan of 
magnetic field lines in both radio sources. The wavelet tadpoles with shorter 
periods (Table ^ show that within the large structure guiding long period 
magnetoacoustic wave trains {P ~ 80s), there are also narrower waveguides 
guiding shorter period waves. 

The complex of radio sources U1-U3 and D1-D3 includes the whole topologi- 
cal structure of the active area contrary to the rather localized EUV sources (Fig- 
ureE]). This is in line with the observed radio sources being usually larger that the 
real sources due to scattering of radio waves {e.g. |Benz, Battaglia, and Vilmer[ 

mm . 

Considering the above presented estimations we can conclude that they sup- 
port the scenario (Figure IH]) of magnetoacoustic waves (wavelet tadpoles) in a fan 
structure above the coronal magnetic null point. 
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